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In  connection  with  design  studies  of  a  high-density-gas-enclosed  Z- 
pinch  fusion  experiment,  Robson  (personal  communication,  1984)  has  noted  the 
advantages  of  surrounding  the  discharge  volume  with  a  water  vortex.  The 
latter  serves  as  a  first  wall,  radiation  barrier,  and  heat  exhange  medium, 
and  scavenges  reaction  products  and  unbumed  fuel.  The  present  note  is 
written  in  order  to  show  how  such  a  vortex  can  be  set  up  inside  an 
appropriately  shaped  convergent  axlsymmetric  duct,  or  infundibulum. 

We  consider  steady  incompressible  flows  in  a  geometry  indicated 
schematically  in  Fig.  1.  A  jet  of  water  is  directed  tangentially  and 
slightly  downward  against  the  surface  of  the  duct  at  the  input  point,  where 
the  radius  is  Rq.  The  duct  axis  is  assumed  to  be  in  the  vertical 
direction.  If  the  inlet  flow  speed  is  v^  and  the  volume  flow  rate  is  V, 
then  the  initial  thickness  A  of  the  layer  of  water  flowing  around  the  duct 
is  given  by 

V  -  2*  RqVqA  sinfl0,  (1) 

where  9q  Is  the  angle  between  the  direction  of  the  input  flow  and  the 
tangent  to  the  circle  of  radius  Rq  in  the  horizontal  plane.  If  A  «  Rq, 
then  all  the  water  has  the  same  values  of  all  three  components  of  the 
velocity  vector  at  the  input  point.  For  the  time  being  we  assume  that  the 

flow  is  laminar. 
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We  can  treat  steady  flows  in  this  configuration  using  either  an 


Eulerian  or  Lagrangian  representation.  Since  the  flow  progresses 
monotonically  downward,  either  distance  traveled  along  a  streamline,  its 
projection  on  the  r  -  z  plane,  or  the  projection  of  the  latter  onto  the 
negative  z  axis  (i.e.,  downward  displacement)  can  serve  as  a  timelike 
variable  parametrizing  the  motion  of  a  fluid  element.  Neighboring  elements, 
however,  although  they  remain  contiguous,  do  not  travel  at  the  same  speed. 

As  a  result,  any  small  volume  of  fluid  undergoes  shear  in  the  course  of  its 
descent  through  the  duct,  so  that  a  Lagrangian  description  based  on  the  use 
of  dlsplacment  Instead  of  time  is  not  very  convenient.  We  therefore  elect 
to  work  in  Eulerian  coordinates. 

The  Eulerian  equations  of  steady  incompressible  axisymmetric  laminar 
flow  in  cylindrical  coordinates  are 


where  we  have  included  the  acceleration  due  to  gravity  in  Eq.  (5). 

To  solve  Eqs.  (2)-(5),  it  is  useful  to  Introduce  curvilinear 
coordinates.  A  streamline  coordinate  a  is  defined  by 


(6) 


da2  ■  dr2  +  dz  , 


where 


dr  u 
Hz’  ”  w  * 


The  unit  vector  in  the  direction  of  increasing  o  is  thus 


-  u-1  (e  u  +  e  w) , 

«M»r  w**«z 


(7) 


(8) 


where 


u  -  <u2  +  w2)1/2.  (9) 

The  coordinate  transverse  to  o  is  t;  it  increases  in  the  direction  of  the  - 
unit  vector 


,-l 


(ew  -  e  u). 


— * 
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(10) 


If  a  -  0  at  the  inlet  position,  then  for  fixed  t,  a  is  the  length  of  a 
streamline  projected  on  the  r-z  plane.  Likewise,  we  define  t  ■  0  at  the 
duct  wall,  so  that  for  fixed  a,  t  measures  the  distance  by  which  a 
streamline  was  separated  from  the  wall  at  the  duct  inlet,  0  <  t  <  d. 

In  terms  of  o  and  T  we  can  express  the  radial  coordinate  of  a  point  on 
streamline  t  by 
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Substituting 


r(o,T)  -  /q da'  u(o',t)/o(0',t). 
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u  1  (w  TC  '  u  TF} 


in  Eq.  (2)  we  find 


9u  +  „  (ilsi  +  |i)  -  0, 
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(12) 


(13) 


(14) 


where  $  -  tan-1(u/w)  is  the  angle  a  streamline  makes  with  the  vertical 
direction,  while  Eqs.  (3)-(5)  become 


3u 

3a 


(ut£  + wf£); 


(15) 


3v  uv 
¥a  “  "  ru ’ 


(16) 


3w 
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(17) 


For  general  inlet  conditions  and  duct  geometry,  Eqs.  ( 14)-( 17)  can  only  be 
Integrated  numerically.  We  can  however,  learn  something  about  the  behavior 
to  be  expected  in  the  general  case  by  considering  a  class  of  flow  fields 
obtained  by  making  some  simple  ansatzes  regarding  the  properties  of  the  flow 
and  then  choosing  a  duct  shape  which  conforms  with  the  solutions.  To  do 
this  we  return  to  the  cylindrical  form  of  the  system,  Eqs.  (2)-(5),  and 
write 


4 


1  3w  .  D. 
T  r  IF  +  r* 


v  ■  C/r; 


w  ■  w(z); 


P  *  P  [A(r)  +  B(z)  j , 


C  and  D  constant,  then  Eqs.  (2)  and  (4)  are  satisfied  identically,  and  we 
are  left  with 


dw  .  .  dB  n 

"  J?  +  *  +  K  ■  °' 


For  Eq.  (22)  to  have  a  solution  it  must  separate  into  two  equations,  one  in 
r  and  one  in  z.  Thus 


1  ,dw,2  1  d2w 

T  W  -  7  w  IF  ’  c<m,t 


1  rdA  C2  +  D2, 

r  ld?  ‘  ~~P - 1  ’  con8t* 


The  solution  of  (24)  is 


w(z)  -  wQ  -  a  (h  -  z)2. 


-  •  ~  «.•  •  •  •  •  -  -  ••  - .  .  • .  >  .  • 


5 


Wq,  a,  and  h  constant,  whence 


I  Ai  -  I  „ 

4  'dz^  2  dz2 


a  w„ 
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and 


A(r) 


K  -  I  (a  wQr2 


.  C2  +  D2* 

+  — rz - ), 


K  constant.  Finally,  Eq.  (23)  yields 


B(z)  -  -  gz  -  w2  -  - 


gz  ~ 


[«0  - 


a  (h  -  z)*] 


Z 1 2 


Equation  (7),  defining  the  projected  streamlines,  yields 


dr  -ct  r  (h-z)  +  D/r 

Hz"  ”  w  -  o  (h  -  z)2  * 
o 


whose  solution  is 


r0  <w(f  a  h2^  +  2Dz 

Wq  -  ct  (h-z)2  ’ 


Tq  constant.  Since  the  duct  surface  must  be  a  streamline,  which  we  can 
specify  by  means  of  Rq,  its  maximum  radius  at  z  -  0,  its  radius  as  a 
function  of  vertical  position  for  z  <  0  is  given  by 


R(z)  -  ( 


-  o  h2)  +  2Dz 
■  a  (h  -  z)^  ^ 


1/2 


Some  remarks  are  in  order  here  regarding  the  limitations  of  this 


solution.  From  Eq.  (26)  we  see  that  v(0)  <  0  implies  Wq  <  a  h2,  while 
w(z)  <  0  as  z  -*  implies  a  >  0.  To  ensure  that  R(z)  >  0  for  all  z,  we 
infer  from  Eq.  (32)  that  D  >  0  must  hold.  In  order  that  the  streamlines 
have  positive  slope  at  z  ■  0,  Eq.  (30)  Implies  that  D  <  a  R^h,  so  that 
h  >  0.  Hence  the  streamlines  have  positive  slope  for  all  z.  From  Eqs. 

(21)  and  (29)  it  is  clear  that  when  -z  is  sufficiently  large  the  pressure 
becomes  negative,  which  imposes  a  constraint  on  the  length  of  the  system. 

And  finally,  we  have  no  way  to  ensure  that  the  pressure  on  the  inner  surface 
of  the  vortex  matches  that  of  the  gas  inside.  Generally  in  experiments  this 
pressure  will  be  roughly  independent  of  z,  which  cannot  be  said  of  the 
present  solution. 

To  make  this  solution  more  meaningful,  we  assign  some  more-or-less 
arbitrary  values  to  the  constants  appearing  in  Eqs.  ( 18— ( 32)  and  examine  the 
consequences : 

Table  1 . 


Parameter 

Value 

Units 

C 

10** 

cm2  s"1 

D 

10 

cm2  s"1 

A 

1 

cm 

h 

102 

cm 

K 

107 

dyne  cm' 

Rn 

10 

cm 

We  give  a  and  Wq  two  different  sets  of  values:  (1)  a  ■  102  cm"1  s"1,  w^  * 
-102  cm  s"1,  and  (2)  a  ■  5  10"2  cm"1  s"1,  Wq  *  102  cm  s"1 .  In  Table  2  we 
show  the  resulting  values  of  the  three  velocity  components  at  z  ■  0;  the 
inlet  flow  speed  vfl;  the  pitch  angle  9-,  given  by 


9 


(33) 


the  duct  angle  ■ 

the  flow  rate  V 

where  we  set  p  < 
(see  Fig.  1)  at 
duct  surface  at 


sin  0n  *  [u(0)2  +  w(0)2]1/2 


|q,  obtained  from  (30)  using 

tan  $0  -  dR/dz  Iq*,  (34) 

given  by  Eq.  (1);  the  corresponding  power  W  *  Vpv2^ 

1.0  g  cm-3;  the  inner  and  outer  radii  R1  and  Rj  of  the  flow 
z  -  -  100  cm  and  z  ■  200  cm;  and  the  pressure  p2  (z)  on  the 
z  ■  0,  z^  *  -100  cm,  and  z2  *  -200  cm. 
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9 


8 
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It  Is  evident  that  for  these  choices  the  dependence  on  gravity  and  the 
radial  dependence  are  both  weak;  moreover,  the  layer  of  water  becomes 


1 

thinner. 

not  thicker,  as  -z  increases,  owing 

to  the 

increase  in  the  vertical 

• 

w 

velocity 

component . 

l' 

b 

„  • 

1 

Table  2. 

■  *  \ 

fl 

Case  (1) 

Case  (2) 

• 

1 

u(0) 

-9,0  cm  s"1 

-24  cm  s-1 

v(0) 

103  cm  s-1 

103  cm  s*1 

1 

w(0) 

-200  cm  s'1 

-400  cm  s-1 

• 

i 

vo 

1.02xl03  cm  s_1 

1.08xl03  cm  s'1 

I 

j 

1 

9o 

11.3  deg 

21.8  deg 

• 

'i 

i 

J 

9 

*0 

2.6  deg 

7.1  deg 

3 

.  i 

I 

• 

V 

12.6*l03cm3  s-1 

25.2xl03cm3s"1 

■ 

• 

w 

655  W 

1.46  kW 

j 

R^Zj) 

6.0  cm 

4.2  cm 

» 

i 

.  ■  -j 

.  j 

R2(zx) 

6.6  cm 

4.7  cm 

.  .  < 
'  .  i 

i 

Rj  (z2) 

4,5  cm 

2.9  cm 

» 

— 

4 

R2(z2> 

4,9  cm 

3.1  cm 

'  r 

•- 

P2(°) 

9,5x10®  dyne  cm' 

-2 

9.5x10®  dyne  cm-2 

.  > 

• 

8,9x10®  dyne  cm- 

.-2 

7.7x10®  dyne  cm-2 

• 

■ 

P2(z2) 

7,4x10®  dyne  cm' 

-2 

5.0x10®  dyne  cm-2 

■ 

» 

•- 

i 

? 

9 

• . 

L 

V  / 

Now  let  us  consider  the  effect  of  restoring  the  viscous  damping  terms 


omitted  from  Eqs.  (3)-(5).  We  assume  that  viscous  processes  are  important 
only  in  a  thin  boundary  layer  next  to  the  surface  of  the  duct,  and  initially 
we  assume  that  the  flow  in  the  boundary  layer  is  laminar,  so  that  the 
kinematic  or  molecular  viscosity  v  can  be  used.  For  flow  in  pipes  this  is  a 
good  assumption  provided  that  perturbations  in  the  flow  at  the  entrance  to 
the  pipe  are  kept  as  small  as  possible  (see,  e.g.,  L.D.  Landau  and  E.M. 
Lifshitz,  Fluid  Mechanics,  Pergamon,  London,  1959,  §29). 

The  thickness  of  a  laminar  boundary  layer  is  given  (Landau  and 
Lifshitz,  op.  cit.,  §39)  by 

5  -  Gt(v  S/v)1/2,  (35) 


where  G1  is  a  geometrical  constant  of  order  unity  and  S  ■  -z  is  the  distance 
traveled  along  the  boundary.  [There  is  some  question  in  my  mind  whether  S 
should  not  be  the  total  length  of  a  (helical)  streamline.]  Taking  Gj  *  1, 

S  »  100  cm,  v  *  103  cm  s”1 ,  and  v  *  10-2  poise,  we  find  that  5  *•  3  10~2  cm. 
The  frictional  force  per  unit  area  of  the  surface  of  the  duct  is  given  by 

a  -  G2p(v  v3/S)l/2.  (36) 


The  power  needed  to  overcome  frictional  losses  in  the  boundary  layer  is 

-  2ir/Ra  vdz 


where  the  bar  represents  an  average. 


2x  vRo  S, 

Taking  the  geometrical  factor  G2 


(37) 

-  1 


and  setting  R  *  5  cm,  we  find  *  100  watts. 

The  boundary  layer  may,  however,  become  turbulent  as  a  result  of  the 
introduction  of  disturbances  at  the  inlet,  through  the  action  of  duct 
surface  irregularities,  or  because  a  critical  Reynolds  number  is  exceeded 


(Landau  and  Lifshitz,  op.  cit.,  $42-44).  If  that  happens,  the  boundary 
layer  thickness  is  given  by 

5  ~  v*  S/v,  (38) 

where 

v  -  G3v^log(v*S/v),  (39) 

Since  vA  varies  slowly  as  a  function  of  S,  5  now  increases  almost  linearly 
with  S,  in  contrast  with  (35).  The  frictional  force  per  unit  area  in  the 
boundary  layer  is  given  by 

o  ■  pv£.  (40) 

Evidently  5  grows  rapidly,  approaching  the  thickness  of  the  whole  vortex  in 
a  few  centimeters  or  less.  When  this  happens,  the  entire  flow  goes 
turbulent  and  boundary  layer  theory  becomes  Inapplicable.  It  is  obvious  in 
any  event  that  viscous  losses  must  be  ouch  larger  in  the  turbulent  case  than 
for  laminar  flow. 


11 


Figure  1 

Schematic  cross  section  of  water  vortex  guided  by  infundibular  duct 
Initial  (z  *  0)  radius  of  layer  is  Rq,  thickness  is  A.  At  some  finite 
z  <  0,  inner  radius  is  R1  and  outer  radius  is  R2.  At  any  point  on  a 
streamline,  n  ■  e_  and  t  -  e„  represent  unit  normal  and  unit  tangent, 

M«w  «w*T  *0*0 

respectively. 


END 

FILMED 

6-85 


DTIC 


